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Abstract 

In  this  paper  the  ajcisymmetric  crack  problem  for  a  nonhomogeneous 
medium  is  considered.  It  is  assumed  that  the  shear  modulus  is  a  function  of  z 
approximated  by  /x  =  This  is  a  simple  simulation  of  materieils  and 

interfaciaJ  zones  with  intentionally  or  naturally  graded  properties.  The  problem  is 
a  mixed  mode  problem  and  is  formulated  in  terms  of  a  pair  of  singuleir  integral 
equations.  With  fracture  mechanics  applications  in  mind,  the  main  results  given 
are  the  stress  intensity  factors  as  a  functions  of  the  nonhomogeneity  parameter  a 
for  various  loading  conditions.  Also  given  are  some  sample  results  showing  the 
crack  opening  displacements. 
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Introduction 


In  solid  mechanics  many  of  the  engineering  materials  such  a^  composites 
and  large  variety  of  bonded  materials  and  structural  components  are  generally 
modeled  as  nonhomogeneous  continua.  The  thermomechcinical  paxcimeters  of  these 
materials  are  usually  assumed  to  be  discontinuous  functions,  very  often  piece-wise 
constants,  in  space  variables.  However,  there  are  some  important  applications  in 
which  the  spatial  variation  in  material  properties  is  continuous.  As  examples  for 
such  materials  we  may  mention  the  following:  (a)  Geophysical  materials  with 
naturcdly  graded  compositions  such  a.s  shale-sandstone.  In  studying,  for  example, 
hydraulic  fracturing  of  such  a  medium  it  would  be  necessary  to  take  into  account 
the  mechanical  property  gradient  of  the  material,  (b)  Interfacial  regions  with 
graded  properties  in  diffusion  bonded  materials,  in  plasma  spray  coating,  and  in 
ion  plating  (Batakis  and  Vogan,  1985,  Houck,  1987).  (c)  Materitd  nonhomogeneity 
resulting  from  temperature  gradients  in  the  solid  in  which  thermomechanicad 
constants  of  the  medium  are  significantly  dependent  on  temperature.  And  finally 
(d)  the  tcdlored  materials  with  predetermined  continuously  varying  volume 
fractions. 

Some  of  the  important  applications  in  the  last  category  of  materials  may 
be  found  in  ceramic  coating  of  metad  substrates  and  in  metal-ceramic  composites 
with  graded  properties.  It  is  known  that  certain  strength  related  properties  of 
ceramic  coatings  can  be  improved  by  layering  the  interfaciad  zone  going  from 
metal-rich  to  ceramic-rich  compositions.  For  example,  in  joining  tungsten  to 
zirconia  by  introducing  four  intermediate  layers  that  contain  80/20,  60/40,  40/60, 
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and  20/80  percent  W/Zr02,  respectively,  it  has  been  shown  that  the  peak  vedue  of 
the  residual  stress  is  reduced  by  approximately  six  folds  compared  to  that 
resulting  from  direct  tungsten/zirconia  bonding  (Hirano,  et  al,  1988,  Hirano  and 
Yamada,  1988). 

The  Fatigue  and  fracture  characterization  of  materials  and  the  related 
analysis  require  the  solution  of  certain  standard  crack  problems.  With  the 
exception  of  the  problem  of  “torsion”  of  an  infinite  medium  containing  a  penny¬ 
shaped  crack  (Kassir,1972),  the  existing  solucions  of  the  crack  problems  in 
nonhomogeneous  materials  with  continuously  varying  properties  have  been 
obtained  under  the  assumption  of  pleine  strain,  generalized  plane  stress,  or  anti¬ 
plane  shear  loading  (see,  for  example,  Dhalival  cind  Sing,  1978;  Gerasoulis  and 
Srivastav,  1980;  Delale  and  Erdogan,  1983,  1988a,  1988b;  Erdogan,  1985;  Konda 
and  Erdogan,  1989).  One  of  the  practical  problems  in  this  area  is  that  of  an 
internal  circular  crack  in  a  nonhomogeneous  solid  subjected  to  “tension.”  The 
problem  is  always  a  mixed  mode  problem.  In  this  paper  we  consider  the  simplest 
of  these  problems,  namely  that  of  an  2uxisymmetric  crack  in  an  infinite  solid  with 
elastic  properties  varying  in  the  cixial  direction  only  (Fig.l).  The  corresponding 
plane  strain  mixed  mode  problem  wcis  recently  considered  by  Konda  and  Erdogan 
(1992). 

Formulation  of  the  Problem 

In  the  crack  problem  shown  in  Fig.l  for  a  large  elastic  solid  it  is  assumed 
that  the  external  loads  as  well  as  the  geometry  are  axisymmetric.  Therefore,  one 
may  easily  separate  the  torsion  component  of  the  problem  in  which  Ug=u(r,  z), 
(Tg,.,  and  (Tgj  are  the  only  nonzero  displacement  and  stress  components.  Similarly, 


by  solving  the  remcdning  problem  under  the  given  external  loads  in  the  absence  of 
a  crack  and  by  using  superposition,  the  problem  may  be  reduced  to  that  of  a 
nonhomogeneous  solid  containing  a  penny-shaped  crack  in  which 


=  pi(r),  a-^,(r,0)  =  P2{r),  0<r<a  (la,b) 

are  the  only  external  loads.  Furthermore,  if  the  crack  radius  a  is  small  in 
comparison  with  the  nearest  distance  from  the  crack  to  the  boundary  of  the  solid, 
then  in  the  perturbation  problem  the  medium  can  be  assumed  to  be  infinite. 
Needles  to  say,  from  the  viewpoint  of  fracture  mechcinics  the  perturbation 
problem  would  contain  all  the  relevcint  information  such  as  the  stress  intensity 
factors  and  the  crack  opening  displacements.  If  we  now  assume  that  locally  the 
elastic  properties  of  the  medium  may  be  approximated  by 


/x(j/)  =  u  =  constant,  (2) 

for  the  axisymmetric  problem  under  consideration  by  using  the  kinematic 
relations 

_  du  .  _u  .  o,  j.  ^ 

~  ^ee  —  ri  ‘■zz  — 


'rr  -  ‘'W  -  ri  -  g^ 

and  the  Hooke’s  law 

<^ij  =  ^ 

the  equilibrium  equations 


^  ^  -0, 

^  +  %‘+k.  =  0, 


may  be  expressed  as  follows; 

(«  +  l){^  -  -  1 


(3) 

(4) 


(5a,b) 


dr2  +  rg^  ^2“  + 


1^}  + 
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(«  +  l)l 


u 


drdz 


+ 


Idu 

^dz 


d^u 


^  \dz 


a^} 


^w\ 

dr'^  / 

+  («  +  l)a|^  =  0. 

(6a, b) 


where  u  and  w  are,  respectively  the  r  and  z  components  of  the  displacement 
vector,  fi  is  the  shear  modulus,  v  the  Poisson’s  ratio,  A//i  =  2i//(l  —  2^'),  and 
K  =  3  —  Au.  In  addition  to  the  boundary  conditions  given  by  (1)  and  the  regularity 
conditions  at  infinity  (requiring  that  u  and  w  must  vanish  as  +  2^-»oo),  the 
crack  problem  must  be  solved  under  the  following  conditions: 


+  <T„{r, +0)  =  cr,.,(r, -0),  0  <  r  <  oo,  (7a,b) 

u(r, +  0)  =  u(r, -0),  u;(r, +0)  =  uj(r, -0),  a  <  r  <  oo.  (8a, b) 


We  now  assume  the  solution  of  (6)  in  the  form 

/CO 

^  F{2,p)p^\{rp)dp, 

/oo 

Q  G{z,p)p3Q{rp)dp,  (9a, b) 

where  Jq  and  are  the  Bessel  functions  of  the  first  kind.  Substituting  from  (9) 
into  (6),  defining  D  =  djdz,  and  inverting  the  related  Hankel  transforms  we  find 

{{k-1]D^  +  q(k-1)D  -  («  +  1)/i2}F  -  {2pD  +  a(«-l)}G  =0, 

{2pD  +  oc{Z  —  k)p]F  +  ■|(/c  +  1)Z)^  +  q(k  +  1)Z?  —  (k  — l)p^|G  =  0, 

(10a, b) 

where  the  following  relationships  have  been  used; 

r)‘^i(^/’)  =  P-^oirp). 
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^Jo(r/j)  =  ^Ji(rp)  -  /)Vo(rp).  (lla-c) 

The  solution  of  the  system  of  the  differential  equations(lO)  is  found  to  be 

k=l 

G{2,p)=  ak{p)Ak{p)e^>‘~,  (12a, b) 

k  =  l 

where  the  functions  A/,,  {k  =  are  arbitrary  unknowns,  are  the  roots 

of  the  following  characteristic  equation 

m*  +  2QTn^  +  {ct^  —Ip^'ya^  —  2p^am  +  ^  a^p^  +  P*  =  0»  (13) 


and  the  coefficients  are  given  by 

2pmk  +  pa{S~K) 


(«:  +  l)m|  +  a(«:  +  l)mfc  —  (/c  — 1)/j 
The  characteristic  equation  (13)  may  easily  be  expressed  as 

(m^  +  am  -  p'^J  +  =  0 

from  which  it  follows  that 


=  (14) 


(15) 


m,=m,=  -  I  +  l|a=  +  +.4,|I^ap}^ 


m,  =  m,  =  -  I  -  l(a"  + 


(16a, b) 


After  solving  for  the  expressions  for  coefficients  given  by  (14)  may  be 

simplified  as  follows; 

„  _  _ 2mfe  a(3  -  «) _  /u-a 

“  o  _L  ■  FTt  w-1  I  \  ’  (^ 

2p  4-  20^(3  -  k)(1  +  K) 

03  =-0-1,  04  ='0*2.  (17a-c) 

By  observing  that  3?(mi,m3)  >  0  and  3f(m2,m4)  <  0,  to  satisfy  the 
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regularity  conditions  at  2  =  ±  00  in  the  solution  given  by  (12)  we  must  delete  the 
terms  involving  and  ^3  for  2  >  0  and  and  for  2  <  0.  Thus,  the  problem 
may  be  considered  cis  that  of  two  half  spaces  involving  two  unknowns,  and  A4  , 
for  the  upper  half  space  and  two  ,  A-^  cind  .43  ,  for  the  lower  half.  Of  these  four 
unknowns  two  may  be  eliminated  by  using  the  stress  continuity  conditions  (7)  . 
Thus,  by  eliminating  /Ij  and  ^43  from  (12)  and  (7)  it  can  be  shown  that 


F{r,p) 


A^e^^^  +  2  >  0, 

(AjAj  +  A2A4)e^^^  4-  (A3A2  +  ,  z  <  0 


(18) 


G(r,p) 


a,A^e^^^  +  a^A^e^^\  2  >  0, 

+  A2A4)e^^~  +  a3(A3,42  +  z  <  0 


(19) 


where  the  expressions  for  the  coefficients  Aj,..,A4  are  given  in  Appendix  A.  The 
unknowns  A2  and  A4  may  now  be  determined  from  the  mixed  bormdary 
conditions  (1)  and  (8). 


The  Integral  Equations 

To  reduce  the  mixed  boundary  conditions  (1)  and  (8)  to  a  system  of 
integral  equations  we  first  define  the  following  new  unknown  functions 

5i(0  +0)  -  u;(r, -0)}, 

52(r)  =  +0)  -  r«(^-0)}.  (20a,b) 

From  (9),  (18),  (19)  and  (20)  it  then  follows  that 

fOO 

9i{r)  =  J  ^  (  -  p{a2  +  G2^2  ~  p{a^  +  G^)A^)pJi{rp)dp, 

roo 

92{r)=  \  {p{1-E2)A2  +  p{l-E4)Ai)pJo{rp)dp,  (21a.b) 

where  the  functions  £?2,“£'4,  G2  and  G4  are  defined  in  Appendix  A.  Equations  (21) 
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determine  A2  and  ^4  in  terms  of  gi  and  g^. 

By  substituting  now  from  (18)  2md  (19)  into  (9)  and  by  using  (3)  and  (4), 
the  boundaxy  conditions  (1)  may  be  expressed  cis 

Iq  Fiip)pJo{rp)dp  0<r<a, 

f  00  1 

J  Q  F2{p)pJ i{rp)dp  =  ■^P2{r),  0  <  r  <  a,  (22a, b) 

Flip)  =  +  "4^4, 

F2ip)  =  V2A2  +  U4A4,  (23a, b) 

where  the  functions  rij,  n^,  V2,  and  are  again  defined  in  Appendix  A.  Inverting 
the  transforms  given  by  (21)  and  from  (8)  and  (20)  observing  that  5fc(r)  =  0  for 

r  >  a,  {k  =  1,2),  the  functions  Fi  and  Fj  are  found  to  be 

A(/>)  =  (-c^U'^i(P«)^i('S)  +  di2Jo{ps)92is))sds  , 

A(p)  =  i  -  d2iJiips)9i{s)  +  d22Joips)92is))sds  ,  (24a, b) 

where  the  functions  d,j(p),  {i,j  =  1,2)  are  also  defined  in  Appendix  A. 

In  order  to  avoid  working  with  divergent  kernels  and  to  simplify  the 
analysis  regarding  the  asymptotic  behavior  of  the  kernels,  first  both  sides  of  (22) 
are  integrated  in  r.  By  using  (24)  it  may  then  be  shown  that 

Iq  diirp)dp^ —  diiJiips)9iis)  +  di2JQ{ps)g2is))sds  =  spj(s)ds  +  Ci^ 

-  ( -d2iJi(p5)5i(s)  +  d22Joips)92{s))sds  =  ^(^^^p2is)ds  +  C2y 

(25a,b) 

where  and  C2  are  arbitrary  constants.  Equations  (25)  are  essenticilly  the 
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integral  equations  that  determine  g-^  and  g2  .  For  /)  -»  oo  it  can  be  shown  that  the 
functions  di-[p)  have  the  following  asymptotic  behavior: 

dnip)  =  -2^  +  0(;,-^),  dM  =  +  0(^-=), 

diAp)  =  +  0(p-^),  dM  =  -  +  0(/)-^).  {26a-d) 

Thus,  from  (25)  and  (26)  we  find 


Q  d"  ^ii(^)"S) )5i('S)  "b  hi2{f‘,s)g2<^s)}sds  —  2^  f  ^  ^Pii^)ds 

g  {^21(^1 d"  (■^22(^5'*)  d-  k22{^,s))g2{s)]sds  —  ^  P2{^)dS'¥C^, 

(27a, b) 

|[if(s/r)  -  £:(s/t-)],  s<r, 


where 


^u(r,5)=  Ji(r/))Ji(sp)(ip  =  ||  ^ 

k  T* 


[if(r/s)-£;(r/s)],3>r, 


fOO  ey  I  f 

^^22(^'S)  =  g  do{rp)Jo{sp)dp  =  ^ 


^K(s/r),  s<r, 


iK{r/s),  3>r, 


(28) 


(29) 


f  00 

f  00 

/ii2(r,s)=  -  Di2{p)JQ{rp)J^{sp)dp, 
r  00 

^21(^5)  =  “  J  Q  ^2iiP)Ji{rp)Jo{sp)dp, 


COO 

h22{r,s)  =  ^  D22{p)JQ{rp)Jo{sp)dp, 


Ai(p) 


^uip) 

D2l{p) 


(du(p) 

^12(/^) 

diiM’ 


_  ‘^2l(p) 
^22(00)’ 


(30a-d) 
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(31a-f) 


On(p)  =  (j 


^22ip) 


-1  , 


22(°°)  :  / 

K  —  1  ^  \  _ 2_ 


‘'i.M  =  -Sfqrt. 


where  K  and  E  are  the  complete  elliptic  integrals  of  the  first  and  second  kind, 
respectively.  By  differentiating,  (27)  can  further  be  reduced  to  the  following 
standard  form: 

•2 


1 

TT 


2 

Q  i  ^  ^  ^2_,(r,s)gj(s)d3  =  ^^^P2(r),  (32a,b) 


where  the  Fredholm  kernels  k-  are  given  by 

Mi(r,5)-1  Mi(r,s)-1 

s  —  r  +  s  +  r 


ku{r,s)  =  ^^-^i.~  +  '  jt;:-  ~  + 


Duip)pJoi^pWiisp)dp, 


CO 


ki2{r,s)=  Di2iP)pJo{^p)JQ{sp)dp, 


00 


fc2i(r,s)  =  7r5  I  ^  D2iip)pJi{rp)Ji{sp)dp, 


^22(^5)  = 


—  ~  ^  A/2(r,s)  —  1 


s  —  r 


s  +  r 


00 


D22iP)pJi{^p)do{sp)dp,  (33a-d) 


Mi(r,s)  = 


M2{r,s)  = 


^(?j,  3  >  r. 


3<r, 


^£(r)  -  5>r. 


(34) 


(35) 


Note  that  the  dominant  kernels  of  the  system  of  integral  equations  (32)  axe 
of  the  generalized  Cauchy  type  (Erdogan,  1978).  Thus,  expressing  the  solution  by 


[a  —  s)  ^s'  * 


(36) 
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and  using  the  function-theoretic  method  described,  for  example,  by  Erdogan 
(1978),  the  characteristic  equations  giving  7^  and  ;?jt  may  be  obtained  as  follows: 


cot7r7fc  =  0,  A  =1,2,  (37) 

cos7r/3i=-l,  cos7r;32  =  l-  (38a, b) 

In  (36)  hi  and  /ij  are  unknown  bounded  functions  which  are  nonzero  at  the  end 
points  s  =  a  and  s  =  0.  Equation  (37)  gives  the  expected  results  71  =  1/2, 
72  =  1/2.  On  the  other  hand,  considering  the  physical  constraints  ^i(O)  =0  and 
u(0,  +  0)  =  u(0,-0)  =  0,  it  may  be  seen  that  the  admissible  roots  of  (38)  are 
/3i=  -  1,  ^2  =  0  giving 


;  _%i/2^l(^)>  52(5)  =  — i/2^2(^)- 


(39) 


(a-s)'/^‘'"  '  (a-s)‘ 

From  (20)  and  (39)  it  can  be  sho%vn  that  for  small  values  of  r/a  the  crack  surface 
displacements  would  have  the  following  asymptotic  form: 


u;  +  _  u; -  ~  u'^  -u~  ^  r  a.  (40) 

Also,  from  the  physical  conditions  rn  —  ru  "  =  0  for  r  =  0  and  r  =  a,  it  is  clear 
that  the  unknown  function  ^'2  nmst  satisfy 


J  I  rg,(r)d 


r  =  0. 


(41) 


In  solving  equations  such  as  (32),  the  accuracy  is  very  highly  dependent  on 
the  correct  evaduation  of  the  kernels  fc,j,  (z,  j  =  1,2).  For  this  it  is  necessary  that 
the  asymptotic  behavior  of  for  s  -»  r  be  examined  and  the  weaJc  singularities,  if 
any,  be  separated  .  From  (26b,c),  (31e,f),  and 

KiX)  -*  log(4/>|l  -A^)  for  A  -♦  1,  (42) 

it  can  indeed  be  shown  that  for  s  -♦  r  the  kernels  k,^  have  logarithmic  singulairities 


-11- 


that  may  be  extracted  as 

— 4^-r  ^  -  f(l -log4^)  +  mii(r,s), 

^  =  ^log  1  ^  I  +  r(^  “  +  '^22> 

c  oo 

fci2(r,s)=  -^asH22{r,s)  -  TTsj^(^Di2(P)P  +  ^yoi^P)Joi^p)dp, 

c  oc 

=f^-s^ii(^,5)  +  7rsj^  (^D2iip)p  +^yi{rp)Ji{sp)dp, 

(43a-d) 

where  Hu  and  H22  given  by  (28)  and  (29)  and  ,  by  virtue  of  (42),  axe  also 
seen  to  have  logarithmic  singularities  and  ttiu  and  17122  axe  known  functions  which 
are  bounded  in  the  closed  interval  0  <  (r,s)  <  a. 

The  Solution 

The  integral  equations  may  be  solved  by  extending  the  interval  and  the 
definition  of  the  unknown  functions  and  the  kernels  into  (  —  a,0)  and  by  using  the 
properties  of  the  Chebyshev  polynomials  T„  and  U„  (Erdogaxi,1978).  These 
equations  may  also  be  solved  by  observing  from  (39)  that  the  orthogonal 
polynomials  corresponding  to  the  weight  functions  Wi{y)  =  (1  +  j/)/-Jl  —  y  and 
W2{y)  =  1/^1  y  =  2s/a-l,  axe  Pi~^'^\y)  and  pi  ^'^\y),  respectively,  and 
by  expressing  the  unknown  functions  Pi  and  §2  as  follows: 

yi('S)  =  <^i(y)  = 

■i^-y  Q 

1  00  r  _i  0^ 

g2{s)  =  4>2{y)  =  \y)  • 

The  integral  equations  may  then  be  regularized  by  using  the  properties  of  Jacobi 
polynomials  and  may  be  solved  by  truncating  the  series  in  (44)  and  using  an 
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appropriate  collocation  technique  (Mahajan,1992).  However,  a  somewhat  more 
efficient  technique  may  also  be  developed  by  defining 

1 

9i{as')  =  .  92{<^s')  =  (45a, b) 

and  by  requiring  that  /i(0)  =  0,  /2(0)  =  0.  The  orthogonal  polynomials  related  to 
the  weight  functions  in  (45)  are 


n(-i3)..^  nn+^)cos{n  +  ^)d  r(n  +  ^)^ 
n\<W  -  n!^ 

(_i  1)  /’(n+i)  /"(n  +  i) 


(46a,b) 


t  =  2s'-i  =  cose,  (i  +  t)f„(f)  =  r„(t)  +  r„  +  i(f),  -i<f<i.  (47) 

We  then  express  and  g2  in  terms  of  the  following  infinite  series: 


^^Pi(as')  =  -  1)  , 

^i92{as')  =  -n7r=n  £  ^2nr„(2y -  1). 

By  substituting  from  (48)  into  (32)  and  by  observing  that 


1  fi  n^tui2s'-i)ds’  _  1  71  r„(t)  +  r„^i(t) 
y  0  N  1  -  s'  s' -r’  ^  J  -  \ 


TT 


(48a, b) 


lfl|  s'  ^n(2^'-  [  ^n-li^)  +  Ur,{x)  Ax\  <  h 

I  C?„(a:)  H- +  i(i)  ,  x>l. 


s'-r' 


r'  =  r/a,  x  =  2/  —  1  ,  UJ^x)  = 


<n 


stn(n+  1)9 
sind 


(49) 


,  I  X  I  <  1,  cosd  =  I,  n  =  0,1,..,  (50) 


(x  —  -Jx^  ~  l) 

G„(x)  =  ,  (x  >  1),  u  =  0,1,... 

■1  X  —  1 


(51) 


the  integral  equations  (32)  may  be  reduced  to  the  following  system  of  equations: 
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where 


+  v,(^)  +  =  P,(x), 

0 

oo 

+  S2nl2U„_i(x)  +  K22„{x)]}  = 

0 

(  —  1  <  X  <  1)  (52a, b) 


Pji^)  =  i  =  l>2,  x  =  ^-l=2r'-l, 

K,i„{x]=^  J'^ki^(ar',as')^Y~,t„{2s'  -  l)ds'  ,  (i  =  1,2), 

=  4  f%,(a'’,a3')^f§f^ds’  ,  (i  =  1,2),  (53a-c) 

j  0  y\s'{l  -  s') 

Also,  by  substituting  from  (48b)  into  (41)  and  by  using  the  orthogonality 


conditions 


'it  =  '1/2,  .  =  ;>0, 

0,  i 


it  may  easily  be  shown  that 


^20  d"  ^-^21  —  0- 


Equations  (52)  and  (55)  are  solved  by  truncating  the  series  and  by  using  a  simple 
collocation  technique.  Very  fast  converging  results  are  obtained  if  the  collocation 
points  Xj  are  selected  as  follows: 

Tn{x,)  =  0i  ''j  =  ^(2i  -  !)■  (56) 


From  the  derivation  of  the  integral  equations  (32)  we  observe  that  the 
right  hand  side  of  (32)  represents  o',j(r,0)  and  cr ,,j(r,0)  for  a  <  r  <  oo  as  well  as  for 
0  <  r  <  a.  Thus,  defining  the  modes  I  and  II  stress  intensity  factors  by 
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^2  =  ^^!?j2(r-a)a,,(r,0), 


(57a, b) 


and  by  using  (51),  from  (47),  (48)  and  (32)  it  can  be  shown  that 

oo  oo 

fcj  =  —  k2  =  —  ■^2n- 

0  0 


(53) 


For  a  homogeneous  medium  (i.e.,  for  a  =  0  in  (2))  modes  I  and  II  crack  problems 
axe  uncoupled  and  the  stress  intensity  factors  are  given  by 


u  — _ A- 

“  7r^^a 


^  7'Pi(^) 
0 


dr  ,  ^2  =  - 


^^3/2 


0  ^P2(^) 
0 


dr 


(59) 


The  Results 

The  main  results  of  this  study  are  the  stress  intensity  factors  calculated  for 
various  loading  conditions  as  functions  of  the  nonhomogeneity  constant  a  by 
defining  the  shear  modulus  in  (i{z)  —  fi.Qexp{az).  Table  1  shows  the  six  different 
loading  conditions  used  in  the  calculations.  The  table  also  shows  the  corresponding 
modes  I  and  II  stress  intensity  factors  in  a  homogeneous  medium  containing  a 
penny-shaped  crack  of  radius  a  obtained  from  (59).  For  the  nonhomogeneous 
medium  the  normadized  stress  intensity  factors  calculated  for  u  =  0.3  are  shown  in 
Table  2.  As  in  the  plane  strain  problem,  for  all  loading  conditions  both  and  ^2 
increase  with  increasing  uq,  the  dimensionless  nonhomogeneity  constant. 
Comparing  the  results  given  in  Table  2  with  the  corresponding  plaine  strain  results 
obtained  by  Konda  aind  Erdogan  (1992)  ,  it  may  be  seen  that  ,  aiside  from  the 
factor  2/t  characterizing  the  difference  between  penny-shaped  and  plane  strain 
cracks,  the  increase  in  the  stress  intensity  factors  k-y  and  k^  with  increasing  aa  in 
the  axisymmetric  case  is  not  as  severe  as  in  the  plane  strain  crack.  For  the  case  of 
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pressurized  cracks  (i.e.,  for  cr„(r,0)  =  -  pg,  (7^,(r,0)  =  0  and  crj,y(i,0)  =  -  pg, 
<Tj.y(x,0)  =  0)  the  compajison  of  the  two  sets  of  results  is  shown  in  Table  3. 

Note  that  the  results  given  in  Table  2  may  be  used  to  obtain  the  stress 
intensity  factors  for  arbitrary  crack  surface  tractions  by  superposition  to  the 
extent  that  the  tractions  may  be  approximated  by  the  second  degree  polynomials 
in  r. 

After  determining  and  p2  or  the  coefficients  B^n  ^Jid  shown  in  (48), 
the  crack  opening  displacements  may  be  obtained  form  (20)  as  follows: 

/a 

9\{s)ds 

+\sin2e') 

ra 

r[u(r, +0)  —  u(r, -0)]  =  —  sg2{s)ds 

T 

=  -  °  '2^  ^^\^B2i{\siv29’  ^'^sivAe')  +  ^  Bjn  ^^5in2n^' 

where 

cos^9'  = 

It  should  be  observed  that  u'^  —  u~  vanishes  for  r  -»  0  .  This  may  be  seen  from 

(40)  and  also  from  (61)  and  (45)  by  using  the  condition  (41)  as  follows: 

lim  [u(r, +0)  —  u(r, -0)]  =  —lim  f[  sg2(s)ds 
r-»0  r-¥Q  J  r 

j  r  a 

=  —  lim  /  sgi(s)ds 
r^OdrJ  r 


Figures  2  and  3  show  some  sample  results  for  the  crack  opening 
displacements.  In  Fig. 2  the  external  loads  are  o’„(r,0)  =  —  Po,  a^,{r,0)  =  0  and  the 
dimensionless  displacement  shown  is  defined  by 


u;(r.  +  0) -u;(r,-0)  (!  +  «)__ 


(63) 


In  Fig. 3  it  is  assumed  that  <T^^(r,0)  =  0,  cr„{T,0)  =  —  Qq  and  the  displacement  is 
normalized  as  follows; 


U{r)  = 


u(r,  +  0)  —  u(r,  -0) 


_(!  +  «:) 
2po 


(64) 


From  figures  2  and  3  it  may  be  seen  that  the  influence  of  the  nonhomogeneity 
constant  a  is  much  more  significant  on  the  axial  displacement  than  on  the  radial 
displacement.  For  reference  we  note  that  for  a  =  0  in  pressurized  plant  strain  and 
penny-shaped  cracks  the  crack  opening  displacements  are  respectively  given  by 


Uy{x,  +  0)  -  Uy{x,  ~  0)]  =  uioNjl  , 

uj(r,  +  0)  -  w{r,  - 0)]  =  -rVa^  ,  u^o  =  '  PoO-  (65a-c) 


In  this  study,  largely  to  simplify  the  analysis,  the  Poison’s  ratio  u  is 
assumed  to  be  constant.  In  an  actual  nonhomogeneous  medium  this,  of  course,  is 
not  possible.  The  assumption  can  only  be  justified  if  the  fracture  mechanics 
parameters  of  interest,  in  this  case  the  stress  intensity  factors,  prove  to  be 
relatively  insensitive  to  variations  in  the  Poisson’s  ratio.  To  give  some  idea  about 
the  influence  of  the  variations  in  v  on  the  stress  intensity  factors,  some  additional 
results  are  given  in  Table  4.  The  table  shows  the  normalized  stress  intensity 
factors  for  various  values  of  i/  and  for  a  small  and  a  large  value  of  the 
dimensionless  nonhomogeneity  constant  aa.  In  these  examples  the  external  loads 
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axe 


assumed  to  be  <T„(r,0)  =  Pi(r)  =  (- po,  -pi(r/a),  -  P2(r/a)^)  ,  0  <  r  <  a, 
<T,,(r,0)  =  P2(r)  =  0.  Particularly  considering  the  fact  that  the  Poisson’s  ratio  of 
the  constituent  materials  and  consequently  that  of  the  composite  nonhomogeneous 
medium  is  likely  to  vary  within  a  much  narrower  range  than  shown  in  the  table, 
the  influence  of  v  on  the  stress  intensity  factors  does  not  seem  to  be  very 
significant. 
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Table  1.  Loading  conditions  used  and  the  corresponding 
stress  intensity  factors  for  q=0. 


Pl(0 

-Pq 

-PiS 

-  P2& 

0 

0 

0 

P2(r) 

0 

0 

0 

-Qo 

-9i(5) 

-  92(5)' 

^1 

^PiW 

0 

0 

0 

0 

0 

■ 

0 

|92>f5 
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Table  2.  The  vaxiation  of  stress  intensity  factors  with  a  for 
various  loading  conditions  shown  in  Table  1,  :/  =  0.3 


<7z.(^.0)  =  Pi{r),  <7„(r,0)  =  0 


.6366 

.6378 

0.6413 

.6465 

.6580 

.6608 

0.6695 

0.6841 

0.7115 

0.7741 

0.8435 


.0 

.0106 

.0212 

.0319 

.0425 

.0532 

.0639 

.0801 

.1073 

.1628 

.2202 

.3413 

0.4714 

0.6105 


0.5 

0.5008 

0.5028 

0.5058 

0.5097 

0.5143 

0.5193 

0.5279 

0.5441 

0.5811 

0.6219 

0.7108 

0.8057 

0.9053 


.4244 

.4249 

0.4263 

0.4284 

0.4309 

0.4341 

0.4376 

0.4434 

.4544 

.4795 

0.5073 

0.5676 

0.6319 

0.6994 


=  0,  (T„{r,0)  =  pg{r) 


0.0 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
•  0.0001 
•  0.0004 


.5 

.5001 

.5003 

.5006 

.5012 

.5018 

.5026 

.5039 

.5068 

.5142 

.5233 

.5447 

.5682 

0.5926 


0.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.0002 


.4244 
.4244 
.4246 
.4248 
0.4252 
0.4256 
0.4262 
0.4271 
0.4290 
0.4341 
0.4403 
.4550 
.4713 
0.4882 


0.0319 

0.0428 

0.0647 

0.0872 

0.1339 

0.1833 

0.2355 


0.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-0.0001 


.3750 

.3750 

.3751 

.3753 

.3756 

0.3759 

0.3763 

0.3770 

0.3784 

0.3821 

0.3867 

0.3976 

0.4097 

0.4224 


Tables.  Comparison  of  the  stress  intensity  factors 
for  pressurized  plane  strain  and  penny-shaped 
cracks  in  a  nonhomogeneous  medium  (i/  =  0.3). 


0 

0.1 

0.25  0.5 

1.0 

2.5 

5.0 

Plane  strain  crack 

1 

1.008 

1.036  1.101 

1.258 

1.808 

2.869 

0.0 

0.026 

0.065  0.129 

0.263 

0.697 

1.567 

Penny-shaped  crack 

0.637 

0.638 

0.644  0.661 

0.712 

0.918 

1.326 

0.0 

0.011 

0.026  0.053 

0.107 

0.280 

0.611 

Table4.  The  influence  of  the  Poisson’s  ratio  on  the  stress 
intensity  factors,  loading:  (r„(r,0)  =  0,  cr^,(r,0)  =  {  —  Po, 
-  Piirla),  Piirlaf),  0  <  r  <  a. 
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t/  a 


Fig.2  z~  component  of  the  normalized  crack  opening 
displacement,  VF  =  (u;"^  —w~)/wq,  Wq  =  (1  +  k)Po®/2Po 
for  the  external  loading  (r„(r,0)  =  —  po,  <7„(r,0)  =  0. 
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r/  a 


Fig.3.  r-  component  of  the  normalized  crack  opening 

displacement ,17  =  (u  —«")/«(,,  Uq  = 

the  external  loading  cr„(r,0)  =  0,  cr„(r,0)  =  —  go- 
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Appendix 

Expressions  for  various  functions  that  appear  in  analysis. 

Ai(/j)  =  (njUa  -  n3U2)/Ai  .  X^{p)  =  (n^uj  -  n3i;4)/Ai , 

>^3ip)  =  («iU2  -  «2Vi)/Ai  .  \^{p)  =  (niU4  -  n4Ui)/Ai .  (AI-4) 

^i(p)  =  (3  -  k)p  -  (1  + /c)mjaj  ,  u_,(p)  =  m_,  +  ,  j  = 

Ai(p)  =  niU3  -  713^1.  (A5-7) 

^  +  p(a3  -  Ui)]  +  (1  +  /c)uj(a3m3  -  aimj)}, 

^i(p)  =  +  I'ajKS  -  K)(a3  -  ai)p 

+  (1  +  K)(m3  -  mi)aia3]},  (;  =  2,4).  (A8,9) 

‘iu(p)  =^(1-^4) -^(1-^:2), 

d22{p)  =  K  -  GJ  - (^2  -  Gj),  (AlO-13) 

^2(p)  —  p((^  ~  •^2)(^4  ~  G4)  —  (1  —  ■£'4)(a2  —  G2)).  (A14) 
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